Closed-form analytic solutions for dilogarithmic double integrals 
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^^■j- ". Abstract. This note presents techniques to analytically solve double integrals of the dilog- 

arithmic type which are of great importance in the perturbative treatment of quantum-field 
! theory. In our approach divergent integrals can be calculated similar to their convergent 

$— i ' counterparts after identifying and isolating their singular parts. 
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1. Introduction and physical motivation. In quantum field theory the differential volume 
element for the three-particle phase space in four-dimensional Minkowski space has the form (ne- 



glecting overall normalization and statistical factors) of a product containing the relativistic four- 
■^J- ■ momenta pi of the three particles: 

o 

g; dPS 3 ~ndW (1) 

23 : 3 

Q-i' Usually one includes energy-momentum conservation (q = Y2i=iPi) and the on-shell requirement 

,_SZh ' for the particles (p 2 = m 2 ) as Dirac delta functions with positive-energy solutions so that 

dPS 3 ^f[d 4 p l 5 + (pl-m 2 )5 A (q-Y /Pl ) , (2) 



i=l \ i=l 



where q is the total momentum of the reaction and rrii are the masses of the final particles under 
consideration. This expression is known as the Lorentz invariant three-body phase space and 
its integration over the kinematic S'-matrix yields numerical predictions for the measurable cross 
sections in physical experiments (see e.g. [1] and [2]). 

The general phase space Eq. (2) has 5 degrees of freedom, resulting from 7 restrictions (3 by the 
on-shell requirement and 4 by energy-momentum conservation) on the overall 12 components of the 
four-momenta. Integrating out all dependence of the undetected spatial components of momentum 
P3, leaves only two degrees of freedom to parameterize the phase space. Therefore, a convenient 
choice to describe the phase space are the two Lorentz invariants [3]: 

y = l 5— and z = 1 5—. (3) 
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(a) 



(b) 

Figure 1. Three-particle phase space. 



°AiV? + A 
(c) 



With this parameterization the three-particle phase space simplifies considerably, namely one finally 
obtains 

d PS3 ~ dy dz, 



and the integration limits introduced in [3] are 



y+ = i-VI 
2y 



z± 



1 - y-i£ + A+^±iv(i-y) 2 -£\/(y-A) 2 -A£ 



(4) 

(5a) 
(5b) 



where the masses squared are denoted as £ = Amf/q with m; = mi = to 2 , and A = rn^/q 2 . Hence, 
the cross sections reduce to the following structure 



JdPS 3 ^2fi~Yl J d v dz fi(y^ z )' 



(6) 



where the individual fi(y,z) are all simple rational functions in y and z. We have found that if 
differential and doubly polarized states are considered [4], these double integrals generalize to the 
types (to, n) 

fi(y,z) 



dy dz 



[(i-y) 2 -^] m/2 [(i-^) 2 -er /2 ' 



(7) 



where to, n = 0,1,2,3 independently. The integrals of (7) will generally display two classes of 
divergences. The so-called collinear singularities manifest themselves as poles and/or logarithmic 
singularities in £ = and occur on the boundaries of the triangle in Figure la. They can be 
regularized by assuming £ 7^ 0, which deforms the phase space according to Figure lb. Here, the 
so-called soft divergences still survive at the origin and take the form of a logarithmic singularity 
in A. Assuming A / finally renders all intermediate phase-space integrals singularity-free. It is 
important to notice that the final cross section (6) is a physical quantity and can not contain any 
divergences. Therefore, in the total sum the limits £ — > and A — ► will always yield finite results. 
The technique of regularization is a means of controlling the intermediate divergences caused by 
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the particular decomposition (7). Typical three-body reactions in quantum field theory consider 
two massive particles and one massless particle, assuming £ / and A = 0. 



2. Convergent integrals. For the convergent case, we use as a demonstration of our technique 
the simplest integral of the type (2, 1): 



f dydz = f 

J [(i-y) 2 -£] va-^-e J c 1 



dy 



In 



(8) 



which contains the phase-space boundaries (5) with A = 0. There is no pole at (1 — y) 2 — £ = 
as the individual contributions I + and corresponding to the limits z + and z_, have similar 
asymptotic behavior with opposite signs and cancel each other. However, for a separate calculation 
of these parts, we require a spurious cut < £ < 1: 

/ = /+ + /_ with I±=± J _ d * _ - In (^(1 ~ z±Y - I + z± - l) ■ (9) 



Next, we consider the substitutions 



fc± = i - y ± ^(i-y) 2 -^ 

for I±, respectively, which gives the following integration limits at the singular points: 

k± = y^±e + 0(^) with 9=y/€^e. 
Applying (10) in (9) and using the new variable v = y/1 — £ in the limits, one gets 



(10a) 



(10b) 



/+ =7f J dk 



{ln£-ln(2-A;)} 



+ 



1 



l+v 



1 



VI 
VI- 



dk 



[k-VC k + VZ. 



ln(2 — k), 



I- =—= \ dk 

VI J 

1-v 



{ln£-ln(2-*0}- 



k-Vl k + Vl. 

Integrals (11) can now be solved without further difficulty in terms of real dilogarithms [5] 

Li 2( ,)=-Ll^. 

J P 

o 



(11a) 



(lib) 



(12) 



Especially for more complicated integrals of the general type (m, n), we use repeated cyclic appli- 
cation of dilogarithmic identities (see also [5]) to dramatically simplify the final solutions. The final 
result is: 



/=2Li 2 
- Li 2 



2 + y/l 
1 + v 



2Li 2 



2 + ^/ +iln Vl- r 



2- V! 

1 + v 



Lio 



1-t) \ / 1 — v 

+ Li 2 



2 + Vf 



In 



£ + 21n(2-^)-±lnU 



1 + u 



Lb 



1 + v 



+ i ln^[ - I ln£ + ln(2 - yft) ] - \ ln 2 (2 - yft) + 



(13) 



It shall be noted, that rigorous numerical tests have confirmed the validity of our method. Full 
analytic results of integrals of this type have been applied in various contexts (see e.g Refs. [3,4] 
and work in progress). 



3. Divergent integrals. For some pathological cases the intermediate phase-space integrals will 
diverge because of additional poles in the integrand functions /j in (7). Then, one needs to regularize 
with A / to avoid contact with the origin as shown in Figure lc. A more intricate example is 
/ = 1/z 2 for integral type (2, 1): 



J 



I \(1 



dy dz 



[(1 - y) 2 - f ] y/(T 



(14) 



Because of the phase-space symmetry y <-> z, it is easy to see that this integral takes the form 

y+ 



2^7 Ja 



dy 1 + 



V(i-y) 2 -er 



(15) 



with the complicated phase-space boundaries (5) for A / 0. For \z\ <C 1, we use the expansion 

z-l + JZ 



In - 



Z-l-yft 



ln[-i^SV^ + (f) 



1 + Vl 



V 



(16) 



to split the integral into 

dy dz 



J 



-I — 

v 2 J JIT 



v+ 



dy 



^(l-^-ey 2 J ^{i-y) 2 -ty 2 V 2V? V z-l-Jl 



1 m(- z - 1 + ^)' z 



(17) 



Ji 



where the second integral J 2 is rendered convergent. Therefore, we can safely put A = in this 
term and proceed as outlined in Section 2 (with two spurious cuts, one for the upper and one for 
the lower integration limit). Straightforward computation yields: 



J 2 = Li 2 



1 + v 
l-v 



Li 2 



l-v 
1 + v 



+ 2 



Li 2 



1^ + Li 2 f 2 l ~^ r 
1 + v V 1 + v- Jl 
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Lio 2 



Li 2 



1 + v + VI 



Li 2 



+ 



+ 



ln£-21n(l- VI) - 12 In i> + § In (jz^) -21n2 In (y^j) 

-| ln£ - In (l + v^) + y + I112I ln£ - \ [± In (l + ^ - In 2] In (l + y/t) 



+ 



Av 
I 



-ln£ + 21n (2- VI) +ilnfi±^j -ln2 + l 



-±ln 2 2 + - 

4 T- 3 



(18) 



The splitting procedure (17) works for all known cases using adequate expansions similar to (16). 
Its main objective is to isolate the logarithmic divergences in A and reduce the difficult divergent 
part of the integrals to integrals of lower type. Here, the first integral of (17) is of the lower 
type (1,0): 

dy dz 1 



/■_! 



V(i-y) 2 -ey 2 ' 

By integrating over z and using (5b), it follows that 

Ji= ItA=M z+ - z -) = 4 I I 



(19) 



A 2 V?+A 



dy V (|/-A) 2 -A| 
4y + £ 



(20) 



which again can be solved by adding and subtracting a splitting term to the integrand. The 
appropriate splitting term is y/(y — A) 2 — A£/£. By neglecting finite pieces in A, one obtains 



Ai^S+ A ° 



1 



+00 



(21) 



The integral form (21) can now be solved by standard integration techniques (e.g. letting t = y — A 
in the first integral). Further expansion in A finally produces 



Ji 



-In A* +iln£ + ln2 + ln(l- y^) - 2 In (2 - yfl) - 1 + 0(*) 



(22) 



so that J = Ji + J 2 in (14) is fully determined. 



4. Conclusion. We have outlined a general integration method to obtain closed-form analytic 
solutions for complicated dilogarithmic double integrals that emerge in applications such as pertur- 
bative quantum field theory. Our method combines suitable substitutions, repeated splitting and 
isolation of divergent integral parts, and cyclic usage of dilogarithmic identities. 

In particular, with this approach it was possible to calculate fully analytic predictions for 
reactions in particle physics [3,4], where in the past only numerical estimates were available. It also 
permitted to find simple and highly accurate polynomial approximations of the Schwinger type for 
these processes [4]. 
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The method has been successfully implemented on computer algebra systems to help create 
extensive tables of relevant integrals of type (m, n). They form the basis of numerical C/FORTRAN 
libraries that will be available in the near future. 
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